Co-tunneling current through the two-level quantum dot coupled to magnetic leads: 

A role of exchange interaction 
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The co-tunneling current through a two-level doubly occupied quantum dot weakly coupled to 
ferromagnetic leads is calculated in the Coulomb blockade regime. It is shown that the dependence 
of the differential conductance on applied voltage has a stair-case structure with different sets of 
"stairs" for parallel and anti-parallel configurations of magnetization of the leads. Contributions 
to the current from elastic and inelastic processes are considered distinctly. It is observed that the 
interference part of the co-tunneling current involves terms corresponding to inelastic processes. 
Dependence of the co-tunneling current on the phases of the tunneling amplitudes is studied. 
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I. INTRODUCTION 

It is well-known that the resistivity of a conductor de- 
pends on its magnetic state There are several quan- 
tum effects increasing this dependence, e.g. Kondo effect 
0|, colossal magnetoresistance effect @, giant magne- 
toresistance effect Qj- One can observe giant magne- 
toresistance effect in thin film structures composed of 
alternating ferromagnetic and non-magnetic layers. At 
the room temperature resistance for the parallel config- 
uration of magnetizations in the leads is a few percent 
less than the resistance for the antiparallel configuration. 
Thus spin degrees of freedom can considerably influence 
the magnetoresistance of the system. 

Recently this influence of magnetism on transport has 
attracted a lot of attention in different systems, in par- 
ticular, in quantum dots [5j-[l5| . A quantum dot is the 
simplest system in which one can study interplay be- 
tween spin and charge degrees of freedom. If the size 
of the dot is small enough the characteristic electrostatic 
energy which is necessary to change the number of elec- 
trons on the dot can exceed characteristic temperature 
and applied source-drain voltage (V). Such regime of 
electron transport is referred to as Coulomb blockade. In 
this regime processes of the first order in tunneling am- 
plitudes (sequential tunneling) are not suppressed only 
at specific values of the gate voltage (so-called Coulomb 
peaks) at which the charging levels for electrons on the 
dot become degenerate in the absence of tunneling [lj| . 
Therefore the current outside the Coulomb peaks is de- 
termined by co-tunneling processes. They are of the sec- 
ond order: tunneling of electrons from one lead to an- 
other through a virtual state on the dot [l?} • 

Addional factor which strongly effects current through 
a quantum dot at low temperatures (T) is discreteness 
of the energy spectrum [18j . For example, in the limit 
T » A the conductance due to inelastic co-tunneling 
is proportional to (e 2 /h)(T/E c ) 2 [l?} where A de- 
notes the single-particle spacing, e the electron charge, 
h the Planck's constant, and E c = e 2 /(2C) the elec- 
trostatic energy of the quantum dot with the total ca- 
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pacitance C. In the opposite regime of low tempera- 
tures T -c A the (inelastic co-tunneling) conductance is 
proportional to (e 2 /h)(A 3 /E 2 T) exp(-ZWT), i.e. expo- 
nentially suppressed with temperature [19j. The same 
holds for the conductance due to clastic co-tunneling: 
at high tem per atures T 3> A the conductance is cx 
(e 2 /h)A/E c |17| whereas for T <C A it is proportional 
to (e 2 /h)(A/E c ) 2 [lj|. We mention that at low tempera- 
tures the elastic co-tunneling contribution dominates the 
inelastic one. 

The presence of exchange interaction of electrons on 
the quantum dot results in dependence of the conduc- 
tance on magnetic polarization of the leads. Recently 
this dependence has been extensively studied in the two- 
level quantum dot !20j. In the paper [2lj numerical so- 
lution of the rate equations describing transport through 
the two-level quantum dot was employed. 

Contrary to papers [2(1 HH, our aim is to calculate 
analytically co-tunneling current through the two-level 
quantum dot with exchange interaction coupled to mag- 
netic leads. Then we examine interference and non- 
interference contributions to the co-tunneling current and 
contributions corresponding to elastic and inelastic pro- 
cesses. We obtain that at non-zero exchange energy 
( J > 0) the interference part of the co-tunneling current 
for antiparallel alignment of magnetization in the leads 
contains inelastic terms which correspond to transitions 
between states of a quantum dot with different values 
of the total spin. This fact differs problem under con- 
sideration from the standard one [17|. Also we demon- 
strate that low temperature asymptotics of the current 
are similar to that of in the two-level quantum dot with- 
out exchange interaction []j[ if we substitute |A — 2J| 
for the level spacing A. At J = A/ 2 when the transition 
from the singlet (J < A/2) to triplet (J > A/2) ground 
state occurs, some of inelastic processes become elastic. 
In addition we find that low-temperature current-voltage 
characteristics are non-linear. However, they became lin- 
ear if the spin-triplet gap | A — 2 J| is equal to zero. This 
fact can be used to observe the spin-triplet transition 
point experimentally. We have found that inelastic co- 
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tunneling current has singularities corresponding to tran- 
sitions between energy levels of the quantum dot. 

In general, contributions to the conductance of higher 
orders (4th and higher) in the tunneling amplitudes con- 
tain logarithmic divergences [22| ■ These divergences [23| 
are not relevant if e\V\,T > T K ~ A exp(-A/V|i| 2 ) 
where v and t stand for the density of states in the leads 
and the tunneling amplitude, respectively. Thus the co- 
tunneling processes provide the main contribution to the 
current. 

To observe the interplay of ferromagnetism and dis- 
cretness of energy spectrum in quantum dots one should 
prepare a quantum dot with large level spacing. In ultra- 
small metallic islands (nanoparticles) [24j or in semicon- 
ducting quantum dots (25| . the level spacing is compa- 
rable with the charging energy or even larger. This also 
happens in the case of a molecule attached to metallic 
leads [26j and in two-dimensional electron gas with fer- 
romagnetic leads (27| . 

The paper is organized as follows. We start with the 
formulation of the model describing the system under 
consideration (Sec. |TTJ) . Then, in Sec. IIII1 we derive gen- 
eral expression for the co-tunneling current under several 
approximations which do not change qualitative proper- 
ties of the system. In Section [IV] we present the expres- 
sions for the non-interference and interference contribu- 
tions to the co-tunneling current. Finally, we present 
discussions and conclusions in Sec. fVl 
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FIG. 1: Eigen states of the isolated two-level quantum dot. 



II. THE MODEL 

The quantum dot coupled to the leads is described by 
the following Hamiltonian: 

H = H,+H r + H QD + H T . (1) 

Here Hj = ej^a^aju, j — Z, r, are Hamiltonians of the 

left and right leads correspondingly with djko / a \ka °- c " 
noting annihilation/creation operators in the correspond- 
ing lead, g - spin index. Term 

H T= E E t( il,aAkJ^+h.C. (2) 

describes the tunneling between the quantum dot and 
lead where d aa stands for an annihilation operator of 
electrons on the dot. The quantum dot is modeled by 
the so-called universal Hamiltonian Q 

Hqd = J2 + E C (N - A ) 2 - JS 2 , (3) 

a,cr 

where N = d\ a d ca denotes the operator of the num- 
ber of particles on the dot, S = J2aa> <$ acr cTo„>d aa i the 
operator of the total spin of electrons, No the equlibrium 
number of electrons on the dot which minimizes the elec- 
trostatic energy and is tunable by the gate voltage, and 
J > the exchange energy. 



The main objective of the present paper is analytical 
investigation of interplay of spin and charge degrees of 
freedom in transport in the co-tunneling regime. The 
simplest system revealing such effects is the quantum dot 
with two single-particle levels. Although Hamiltonian Q 
is derived under assumption that the number of levels on 
the dot involved in the transport is large [9|, it is wi dely 
used to model a two-level quantum dot (20|, HI], US l29j ). 
As we demonstrate in Appendix [X] Hamiltonian can 
adequately describe the many-particle spectrum even for 
a quantum dot with only two low-lying single-particle 
levels. 

In what follows we consider two-level quantum dot, i.e. 
e a will take only two values t\ and ti = e\ + A. There 
are 16 (many-particle) eigen states. The energy of a state 
denoted by \NSnM) is given by the following expression 



E NS nM = E c (N-N ) 2 +ne 1 + (N-n)e 2 -JS(S+l). (4) 



Here N denotes the number of electrons on the dot, S the 
total spin of the dot, and n the number of electrons on 
the first level. The eigen states of the isolated quantum 
dot are depicted on Fig. [TJ 
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III. CONTRIBUTION TO THE CURRENT OF 
THE FOURTH ORDER IN t[ j) , 

A. Perturbation theory 

Let us assume that the Coulomb energy is much larger 
than other energy scales, E c ^> T, A, J, and N is close to 
the integer number (regime of a Coulomb valley). Then 
the main contribution to the current at low tempera- 
tures is of the 4th order in a a' E3 • Indeed the 
odd terms vanish. The second order term correspond- 
ing to real transitions changing the charge of the dot is 
exponentially suppressed (~ exp(— E c /T)) for almost all 
gate voltages except those which correspond to Coulomb 
peaks [l6| • The contribution to the current of the fourth 
order describes transitions of an electron from lead to 
lead trough the virtual states on the dot. These transi- 
tions are referred to as co-tunneling of electrons. 

Let us consider the equation for the density matrix in 
the interaction representation: 

p=-i[H T ,p}. (5) 

By solving it perturbatively, we find for the third order 
correction [30l ] 

/t 3 pt 2 ptl 

dt 2 / dh / dt 
-oo J — oo J ■ — oo 

x [H T (t 2 ),[H T (h),[H T (t ),p ]}] (6) 

where po = exp(— /3(if; + H r + Hqd)) stands for the 
cqulibrium density matrix of the unperturbed system 
(J3 = 1/T). The current operator is equal to derivative 
of the charge operator in one of the leads: 

Ij = e4 =iX- lX\ X = 4*;aa<4fc<Aa'. ( 7 ) 

kaaa' 

It is worth to mention that the current is one order higher 
in tunneling amplitudes than the density matrix. The 
fourth order in t j , correction reads 




x (pt\t)H T {t{]H T (t 2 )H T {U)) 

- (E\ T {H)X\t)H T {tr)HT{tJj) 

- {fi T (t 2 )X\t)H T (t x )H T {t 3 )) 
+ {H T {t 3 )H T (t 2 )X\t)H T (t{)) 

- (HT(h)XHt)H T (t 2 )H T (t 3 )) 
+ (Ht (t 3 ) Ht (ti)X^(t) Ht (i 2 ) ) ) 
+ (H T (t2)HT(h)Xl{t)H T (t 3 )) 

- (HT^HT^HTit^xHt))^ (8) 



where Z = Tr(po) is the grand canonical partition func- 
tion and (...) = Tr(. . . po)/Z. 

Let us express all of the averages in Eq. (JS) in terms 
of the exact two-particle correlators for the isolated dot 
(Hqd) and the following Green functions of electrons in 
the leads: 

Here we introduce /? = {fctr}. Each of the eight terms 
in Eq. flS) is evaluated in Appendix (cf. Eqs (|BTj) - (|B7|l ). 
For example, 

(X\t)H T (h)H T (t 2 )H T (t 3 )) 

= (di 1 tda 2 t 1 dl l3t2 d a4 t 3 ) T 1 4 32 G l ^ ij34 (t,t 3 )G r ^ j32 (t 2 ,t 1 ) 

- (di lt da 2 t 1 d'l ;jt2 d ai t 3 ) t 1234 G 1 ^ iI32 (t, ti)G r ^ 0i (t 2 ,t 3 ) 
+ (d^tdi^da^daits) Ti 32A G l p ih (t,t 2 )G r ^ f ^{ti,t 3 ) 

~ (da^dl^dastida^) ^1423^^^ (t, t 3 )G£>j 3 (il , t 2 ) 
+ (da 1 tda 2 t 1 d a3 t 2 da il t 3 ) TL^G^^ (t, h)Gp<p 3 (t 3 , t 2 ) 

- (di lt d a2 t 1 d a3 t 2 dl lit3 ) Tia42G l g* 1 p s (t,t2)Gp<p a (ta,ti), 

(10) 

where (...) = Tr . . . e"^ 13 / Tr e ~ /3ff « D , a k = {a, a}, 

Tijkl = Wi^cy^aAa,' and d «* = d *( f )- We mei> 

tion that Eq. (|10l) involves terms which are proportional 

t0: V iai 42«2*ft« 3 4a4 aild tpiajfcajpsastpiai- If WC 

set, e.g. V = 0, the current should vanish. However, 

the terms proportional to % iai t l p 2a . 2 % 3 a 3 tl p iai remain un- 
changed. Thus such terms give no contribution to the 
current in the fourth order. Therefore we shall omit them 
in what follows. 

Due to the presence of interactions in Hqd, correlators 
of the form (d^dd^d) in Eq. (| 10[) cannot be simplified 
with the help of the Wick theorem. In general absence of 
the Wick theorem leads to a very tedious expression for 
the current. Therefore we introduce some simplifications 
which do not affect qualitative properties of the system 
but allows analytical calculation of correlators (d^dd^d). 



B. Approximations 

We calculate the 4th order correction to the current, 
Eq. (J5J under the following assumptions. 

(i) The leads are made of ferromagnetic metal with 
magnetization along some axis z. Considering 
the exchange interaction to be isotropic we obtain 
that the Green function of electrons in the lead: 
G fcffc 2CTl( r 2 (*>0 = -iteikiaiWaJfcoJO) is Propor- 
tional to A(ki, k 2 ,t, t')4ri,<72 + B(kx,k 2 ,t, t')a% i (J2 , 
i.e. it is diagonal in the spin indices a\,a 2 . Under 
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this assumptions Green functions read (j = I, r) 



x e -^-*')<5( £ _ efciCTi ) 

= -iSk^S^zg&g^tJ), (11) 

= * 1 ^ ffl „ 2 4?C(^')- ( 12 ) 

Here t ka stands for the energy of a single-particle 
excitation in the leads and coefficient Z& arises 
due to renormalization of the spectral density due 
to interactions in the leads. 

(ii) Next we assume that there are no magnetic im- 
purities in the tunneling junctions. Therefore we 
neglect the probability for the spin to flip during 
the tunneling, i.e. 



t 



(13) 



Also since only energies near the Fermi level in the 
leads are essential for the calculation of the current 
we will ignore the dependence of the tunneling am- 
plitudes on energy and introduce (dimcnsionless) 
tunneling conductances 



1 

A ^ 

k 



Z^8{E F -e k *)\t 



kea I 



(14) 



(iii) Finally, we restrict our consideration to the states 
with only two electrons on the dot. If A > 2J 
the spin in the ground state is zero (both electrons 
occupy the lowest energy level). In the opposite 
case of A < 2 J the ground state is ferromagnetic 
(total spin is unity) such that each energy level is 
singly occupied. Therefore it is possible to observe 
drastic effect in the transport through the quantum 
dot due to the change of the spin in the ground 
state. In what follows we assume such value of the 
gate voltage that the only important states of the 
system are the states with two electrons on the dot, 
i.e. N ~ 2. 

In the current |(5J) one can distinguish terms of two 
types. Terms of the first type depend only on the ab- 
solute values of the tunneling amplitudes (cx |i'| 2 |': r | 2 ). 
Contributions of the second type involve also relative 
phases of the tunneling amplitudes. One can refer to 
corrections of the first type as non-interference whereas 
the second type as interference contributions. To single 
out effects associated with dependence of the current on 
the phases of the tunneling amplitudes it is convenient to 
analyse interference and non-interference contributions 
to the current separately. As well-known [l7| . there are 



processes of two types: inelastic processes during which 
the energy of the quantum dot changes and elastic with 
the same energies of the initial and final states of the 
quantum dot. Due to conservation of the energy of the 
whole system we have: e\ + Ei = e 2 + Ef, where t\, e 2 
denote energies of an electron before and after the tun- 
neling event, Ei, Ef energies of the initial and final states 
of the dot. By definition inelastic co-tunneling involves 
change of the energy of the dot: Ei 7^ Ef, (electron-hole 
pair arises). It means that during the inelastic cotunnel- 
ing the state of the quantum dot changes \i) 7^ |/). So 
one can call this process as incoherent [3l[. Elastic co- 
tunneling process can either change the state of the dot 
(\i) 7^ I/) the electron spin flips) or does not \i) = \f). 
One can refer to the latter type of co-tunneling as co- 
herent. In what follows each term in the current will be 
discussed according to the definitions introduced above. 



IV. CO-TUNNELING CURRENT 



A. General expression 



Working out Eq. ((5J) we obtain the following expression 
for the current 



( 4 ) _ _2_ F cfei&g f 
nm ~ Z 2 J (2tt) 2 



(l-n^( £l ))^( e2 )x < ' > (£i,e 2 )+ 



+ (1 - n l F (ei))(l - n F {e 2 )) X <<<{e 1 ,E 2 )+ 
+n F (e 1 ){\~n F {e 2 ))x > ' < {e 1 ,e 2 )+ 

+n F {e 1 )n r F {e 2 )x > ' > {e 1 ,e 2 )\. (15) 



Here n F (e) stands for the Fermi-Dirac distribution of 
electrons in left/right lead. Terms in Eq. (fT5j) propor- 
tional to (1 — n l F (ei))(l — n r F (e 2 )) and n l F (ei)n F (e 2 ) have 
no physical meaning. Therefore x > ' > an d X <,< have to 
vanish; it is in agreement with direct calculations. The 
detailed expressions for x <,> an d X >,< arG cumbersome. 
As an example, we present expressions for parts of x <!> 
and x > ' < which contribute to the non-interference part 
of the co-tunneling current are presented in Appendix [B] 
(cf. Eqs. (|B9)l - ((BT2"|) ). Next, 



Z 2 = e-^ ei+e2 \t 



/3A 



-/3A 



+ 3e 2/3J + 1) (16) 



is the canonical partition function for the Hamiltonian 
Hqd with N = 2 electrons. Analytical expressions for 
non-interference and interference contributions to x's are 
presented below for antiparallel and parallel alignment of 
magnetization in the leads. 



B. Antiparallel configuration 

Antiparallel configuration of magnetization in the leads 
corresponds to the following values of g l ' r : 



ffi.T = ff2,t = s'j 9\,±= 9l,x = 9 r , 
9 l i,i = .924 = fi,t = -92 : t = °- 



(17) 



This choice of the values of g l ' r assumes that spin-down 
electron band of left lead and spin-up electron band of 
right lead are empty. In addition we introduce phases <pi 
and d) r : 



<"lt 2f 



^l\\\^2\\ e 



(18) 




In the anti-parallel case we obtain (<f> = 4> r — 4>i) 



47TAVg^ .-ff( ei+ e 3 -2J) 

El 

2S(E! -e 2 )(l + cos^) 
-8{ei - e 2 + 2J)(1 - coscj)) 
-6(e x - e 2 - A + 2 J) 

-6(5! ~ £ 2 + A + 2 J) 

_ -2p\7 



<5(ei — £2 — 2J)(1 — COS0) 
+e ^(A-2j )j(£i _ £2 + A _ 2 J) 

+e -m+^s{e 1 -e 2 -A-2J)], (19) 



and 



£2) = -Xap (£2,ei)- 



(20) 



Each term in Xap" nas transparent physical interpreta- 
tion; it can be written as ~ exp(—Ei/T)6(ei — £2 + £/ — 
Sj). For example, the term proportional to 6(61 — 62) cor- 
responds to the elastic cotunneling of an electron which 
results in the transition of the quantum dot between the 
state |2110) and [211 — 1>. Similarly, the term in Eq. dT5J) 
proportional to 6(ei — £2 — A + 2 J) describes inelastic co- 
tunneling and corresponds to the transition of the quan- 
tum dot from the state 1 2 1 1 — 1) to 1 2020) . This transition 
can be realized through two virtual states: with three 
electrons (upper path) and one electron (lower path) on 
the dot. Under our assumption of large E c both of them 
provide equal contribution to the current: the energies 
of the virtual states are equal to E c with our accuracy. 
Non-zero value of exchange energy J allows the terms in 
X <!> (proportional to 6(e± — £2 ±2 J)) which are inelastic 
but depends on phase difference (f>. They correspond to 
transitions between the states 1 2 1 10) and 1 21 1 ± 1) (see 
Fig.©. 



FIG. 2: Example of inelastic co-tunneling process that gives 
contribution to the interference part of the current for anti- 
parallel configuration. See text. 



Utilizing Eq. (fTo) we find 

T< & = "^fr e-^ 1+e2 - 2J) [2F(U)(l +cos0) 

+ F(V + 2J)(1 - cos 0) + F(V - A + 2 J) 
+ F(V + A + 2J) + e-' 2 P J F(V - 2J)(1 - cos, 
+ e m-2J) F (y + A-2J) 

+ e -^ A+2J ^F(V — A — 2 J) — (V -» —V) 



_ Mnel) T {el) Mnel) j(el) 

AP pi™ trim T 1 4 P in ' 1 A P 



l AP. 



APA 



APA 



(21) 



where F(e) = e/[exp(e/T) — 1]. The inelastic and clastic 
contributions of non-interference part of the co-tunneling 
current are given as 



(inel) 



2A 2 g l g r 



AP t nin 



-/3(ei+e 2 -2J) 



F(V + 2J) 
A + 2J) 



^Z 2 El 
+ F(V- A + 2J) + F(V 
+ e- 2f3J F(V - 2J) + e /<- A ~ 2J ^F(V + A - 2 J) 
+ e ^- A - 2J ^F(V - A - 2J) - (V -> -V)] , 

(22) 



and 



T (el) 

AP,nin 



4ff ff r A 2 
ttESZ* 



/3(£l+£ 2 -2J) 



V, 



(23) 



respectively. Here we use the following relation: F(V) — 
F(—V) = —V. The interference term of the co-tunneling 
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current is splitted on inelastic and elastic parts as follows 

2A 2 5 y 



jiinel) 
l AP,in 



nZoE? 



-/3(ei+e 2 -2J) 



F(V + 2J) cos< 



+ e- 2/3J F(V- 2J)cos0- (V 



-V) 



and 



I 



W _ 4gVA 2 , (£i+£2 . 



2./) 



Fcosi 



(24) 



(25) 



As we have mentioned above some inelastic transi- 
tions (e.g. from 1 2010) to |211 — 1)) can be implemented 
through two different virtual states. The interference of 
these two processes depends on the phases of tunneling 
amplitudes. Therefore there is the interference term in 
the current which involves inelastic contibutions of the 
form F(±2J + V) (see Eq. flU}). In the regime of low 
temperatures and voltages, |V|,T, |A— 2 J| <C A, J the in- 
clastic terms of the interference contribution (j2"4")l are sup- 
presed due to small exponential factor exp(— 2J/T). In 
the case of J = the contribution I^p^n becomes elastic 

and exactly compensate I^pm- Therefore, at J = the 
co-tunneling current becomes independent of the phase 

<f>- 

In Fig. [3] we present the dependence of the differen- 
tial conductance dl^p/dV on T and V for 4> = ir/2 
and different values of the exchange interaction J. At 
low temperatures, T « J,A the stair-case structure in 
differential conductance appears. It corresponds to the 
step-by-step switching on of the different inelastic pro- 
cesses with increasing voltage. As follows from Eq. (f2Tj). 
one can expect features at \V\ = |A — 2J|,2J, A + 2 J. 
However some of steps have exponentially small height at 
low temperatures and, therefore, are invisible as shown 
in Fig. [3] At J = there is only one stair-case which 
corresponds to change of the energy of the quantum dot 
by AEqd = A. For small values of the exchange energy, 
J < A/2, the feature at \V\ = AE QD = A-2J is visible. 
In the regime A — 2 J <C A, J, the stair-cases correspond- 
ing to processes with AEqq = A + 2J, AEqd = 2J 
and AEqd = A — 2J appear. As expected, the latter 
disappears at J = A/2. All three stair-cases survive 
at J > A/2. The evolution of the differential conduc- 
tance dl^p/dV with increasing J at fixed temperature 
is shown in Fig. @]for <j> = ir/2. At J > A/2 the feature 
corresponding to the inelastic process with AEqd = 2 J 
disappears at = 0. 

In the most interesting regime of vicinity of the singlet- 
triplet transition where A = 2 J — k with |V|,T, \k\ <g; 



(inel) 



A, J the expression for I Ap 
fled: 



can be drastically simpli- 



(inel) 
AP,nin 



2g l g r A 



tt(3 + e- K l T )El 



(V + k) 



1 



oV/T 



1 _ e (V+K)/T 



+ (V-k) 



1 



-V/T 



1 - e (-V+K)/T 



(26) 



At k = the current ([26]) acquires particularly simple 
form: / 



inel) 
AP,nii 



g g r A V / {kE 2 ), since the spin flip pro- 
cess becomes elastic. This is the reason why the conduc- 
tance does not turn into zero when V = 0. For |k| 3> T, 
the curent I^pnin becomes 

T (,inei) 2c K A 2 g l g r M/T 

1 AP,nin ~~ — " e 



where 



'2(|«;|sh^ + y[l-ch^]), 
e l*l/ T (V-|K|sgn(tO), 



1/3, k>0, 

1, K<0. 



\V\ < \k\ 
\V\ > \k\ 



(27) 



(28) 



Expression (|27[) demonstrates exponential suppression 
of conductance at low temperatures due to the spac- 
ing between energy levels (in our case it is spacing be- 
tween triplet and singlet energy levels which is equal to 
|A-2J|). 



C. Parallel configuration 

Parallel configuration of the magnetization in the leads 
corresponds to the following values of g l,r : 



9i,l = 9 2 ,i = 9 i,i = 92,i 



0. 



This choice of the values of g l ' r assumes that spin-down 
electron band of both leads are empty. In addition we 
introduce phases <f>i and <fi r : 



(30) 



In this case we obtain 



Xp > (ei,£2) = - 



2^A 2 5 V 



-/3(ei+e 2 -27), 



El 

2«5(ei -e 2 )(l + cos0) 
+5(ei - e 2 - 2J)(1 - cos0) 
+6{ei -e 2 -A-2J) + S( £l 



2J) 



+e 



-2/3J 



5{S\ — E2 + 2J)(1 — COS( 



+8{e 1 - e 2 - A) + 5{ei - e 2 + A) 



J(A-2J) 



5(ei - ea)(l - cos. 



-S(S! -S2-A + 2J) + 5( Sl - e 2 - A) 



+e -/3(A+2J) < 5( ei _ e2 )(l_ COS 0) 

+S(ei - e 2 + A + 2J) + 8(e 1 - e 2 



A)]}, 
(31) 
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FIG. 3: Dependence of the differential conductance dl^'/dV on V and T for the antiparallel (left) and parallel (right) 
configurations for different values of J. We use gi = 0.1, g r = 0.2, E c /A = 10, cfr = ti"/ 2 and J/A = 0.2, 0.45, 0.5, 0.55, 0.95 from 
top to bottom. 
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FIG. 4: Dependence of the differential conductance dl^/dV 
on V and J for the antiparallel (left panel) and parallel (right 
panel) configurations at (j> = 7r/2. We use gi = 0.1, g r = 0.2. 
E c /A = 10 and T/A = 0.01. 



and 



Xp < (ei,£2) = ~Xp > (£2,£2) 



(32) 



As in the case of the antiparallel alignment of mag- 
netizations, each term in Xp' > nas transparent physical 
interpretation. As compared with \ , Eq. (f3"T1) demon- 
strates that the case of parallel magnetizations allows 
more elastic processes. In addition, transitions in which 
energy of the quantum dot is changed by AEqd = ±A 
are possible. 

Utilizing Eq. (fT5"|) . we find 



(4) 



AVfl r 



-/3(ei+£2-2J) 



+ F{V + 2J){\ -cos( 
+ F(V + A + 2J) 



{2F(V)( 



1 + COS ( 



F(V- A + 2J) 



F(V-2J)(l-cos< 



i^(l^ + A) + F(V — A) 



+ e 0(A-2J) 



F(V)(l -cos( 



+ F(V + A - 2 J) + F(V + A) 

+ e-P {A+2J) \F(V)(l-cos4>) 

+ F(V-A-2J) 

+ F(V — A)j -(V->-T0} 



j-(inel) 



Ael) 



j-(inel) 



,(ei) 



(33) 



P,nin ' P^nin ' P,in ' P,in' 

The inelastic and elastic parts of the non-interference 



contribution to the co-tunneling current are as follows: 



j-(inel) 
l P,nin 



-/3(ei+e 2 -27) 



|f(U + 2J) 



F(y — A + 2 J) + F(V + A + 2 J) 



(34) 



-2/3 J 



F(7 - 2 J) + F(y + A) + F(V - A) 



P /3(A-2J) 



-,S(A+2.7) 



F(V + A - 2 J) + F(V + A) 
-F(V-A-2J) 
F(V — A) -(V->-V)} (35) 



and 



P.nin 



2A 9 l 9 r 0(e 1+ e 2 -2J)( 2 

nZ 2 E* \ Z 



-0(A+2J) 



o/ 3(A-2J) 



V. 



(36) 



The inelastic and elastic terms for the interference part 
of the co-tunneling current are given by 



j-iinel) 
l P,in 



AV5 



ttZoE? 



-f>{ 



- 1+t2 ' 2J ^F{V + 2J) 



-2/3J 



f(v — 2 j) 



(V -> -F)}cos0 



(37) 



and 



r ( e {) 



_2gVA2 

e -/J(A+2,/)j Fcos ^ 



-/3(ei+6 2 -2J) 



2 _ e /3(A-2J) 



(38) 



In the expression (|33j) there arc two types of additional 
terms in comparison with the case of the antiparallel 
configuration. First ones (oc F(V ± A))) correspond to 
the singlet-singlet transitions of the quantum dot (e.g. 
1 2020) — > |2010)) during inelastic co-tunneling including 
a transfer of one electron to the different level. Second 
one (oc F(V)(1 — cos<fr)) are elastic terms due to tran- 
sitions of the quantum dots between the states in which 
one-level is doubly occupied, e.g. |2020). These elastic 
terms lead to dependence of the co-tunneling current for 
parallel alignment of magnetizations on the phase differ- 
ence <f> even at J = 0. 

In Fig. [3] we present the dependence of the differen- 
tial conductance dip /dV on T and V for (f> = n/2 
and different values of the exchange interaction J. At 
low temperatures, T <C J, A the stair-case structure in 
differential conductance appears. It corresponds to the 
step-by-step switching on of the different inelastic pro- 
cesses with increasing voltage. As follows from Eq. (|33p , 
one can expect features at \V\ = |A — 2 J\, 2 J, A, A + 2 J. 
However some of steps have exponentially small height at 
low temperatures and, therefore, are invisible as shown 
in Fig. [3] At J = there is only one stair-case which 



9 



corresponds to change of the energy of the quantum dot 
by AEqjj = A. For small values of the exchange energy, 
J < A/2, additional feature at \V\ = AE QD = A - 2 J 
is visible. In the regime A — 2 J <C A, J, the stair- 
cases corresponding to processes with AEqd = A + 2 J, 
AEqd = 2 J and AEqd — A — 2 J appear. As expected, 
the latter disappears at J = A/2. As in the case of an- 
tiparallel alignment of magnetizations three stair-cases 
at \V\ = 2J - A, 2 J, 2 J + A survive at J > A/2. The 
feature corresponding to AEqd = A becomes invisible 
at J > A/2. The evolution of the differential conduc- 
tance dl^p 1 jdV with increasing J at fixed temperature is 
shown in Fig. [4] for cj> = n/2. At J > A/2 the feature 
corresponding to the inelastic process with AEqd = 2 J 
disappears at = 0. 

In the most interesting regime near the singlet-triplet 
transition A = 2 J — k: \V\,T, \k\ <C A, J the expression 



for Ip^n can be written as 



r(inel) 



A 2 g l g r 



7r£;2(3 + e - 
1 - 



+ {V + k) 



e V/T 



(V-k) 



I _ e (V+ K )/T 



1 - e- v ' T 

1 - e {-V+K)/T 

(39) 



(l/2)I { ™ el) 



We mention that in this regime I Pnm — v-v ")± AP,nin- 
In the case |«|, |V|, T <C A, J additional terms oc F(V ± 
A) are suppresed by the small factor exp(— A/T) and 
therefore they do not contribute to the current. 



V. DISCUSSIONS AND CONCLUSIONS 

In this paper the co-tunneling current through the two- 
level quantum dot coupled to ferromagnetic leads is cal- 
culated analytically. The results have been presented for 
the most interesting case of the quantum dot with two 
electrons and for parallel and antiparallel configurations 
of magnetization in the leads. 

Inelastic co-tunneling current has features correspond- 
ing to transitions between energy levels of the quantum 
dot which could be used to determine structure of en- 
ergy levels on the quantum dot. In both cases of par- 
allel and antiparallel configurations the non-interference 
part of the conductance has a minimum at low temper- 
atures and voltages |V|,T <C A. Width of this mini- 
mum near the transition (A = 2 J) between singlet and 
triplet ground states is defined by the gap between sin- 
glet and triplet two-electron states |A — 2J|. Near the 
transition our results for the differential conductance re- 
sembles the expression derived in the paper [l9| for the 
two-level quantum dot at J = 0, if we substitute the av- 
eraged single-particle spacing by the singlet-triplet gap 
| A — 2J|. In the case of low temperatures T -C |A — 2J\ 
inelastic part of the current is suppressed by the factor 
cxp(— | A — 2J\/T) in comparison with elastic one. How- 
ever at temperatures T ~ |A — 2J| <C A elastic and in- 
elastic parts of the co-tunneling current are of the same 



order of magnitudes. It is due to the fact that the spin flip 
processes become almost elastic near the singlet-triplet 
transition in contrast to the case of J = 0. 

As we mentioned in the Introduction numerical calcu- 
lations of the differential conductance based on a rate 
equation approach were performed in Ref. [2l[. It was 
found that there is a zero-bias peak with the width of 
the order of T in the antiparallel case (in contrast to the 
parallel one). This zero-bias peak has been explained by 
non-equilibrium difference in the occupation probabili- 
ties of states |2111) and |211 — 1). The absence of such 
effects in perturbation theory may be a reason for dis- 
cordance between our findings and results of Ref. [2l| at 
small bias. However, we emphasize that both perturba- 
tion theory and rate equation method produce the same 
results for conductance and for the positions of "stairs" 
in dependence of the differential conductance on voltage. 

In the regime V, T ~ | A - 2 J| < A only singlet-triplet 
transitions are important. It is worthwhile to mention 
that in this regime the inelastic interference contributions 
to the co-tunneling current arc cxponentialy suprcsscd. 
Thus, in this regime expressions (|2T[) and (|39|) for inelas- 
tic part of the co-tunneling current j( me ') are valid for 
a quantum dot with large number of levels. Elastic part 
of the co-tunneling current 1^ in such quantum dots 
is determined by transitions through energy levels in the 
range ~ E c . Therefore, expression for 1^ in quantum 
dots with large number of levels is ~ E c / A times greater 
than 1^ for a two- level quantum dot (see Eqs. (|2"3")l . 
(|25|) . (|36|) and ([38])). It is this way the latter matches 
with the result for the clastic co-tunneling in a multi- 
level quantum dot. 

In quantum dots with large number of levels tran- 
sitions between ground states with 5 and 5 + 1 are 
possible at J = J s = A(25 + l)/(25 + 2) [|. Our 
results indicate that in vicinity of such transitions at 
\V\,T < | J - J$\ < J, A the inelastic part of the co- 
tunneling current will be suppressed. However, exactly 
at the transition (J = Jg) /( me/ ) become linear in V 
and independent of T. Therefore, at the transition point 
the increase of conductance should occur. One can utilize 
this fact to experimentaly observe the transition between 
ground states with 5 and 5 + 1 in multi-level quantum 
dots. 

It is worthwhile to emphasize that the interference part 
of the co-tunneling current (|2^|) and (|3"T|) involves terms 
corresponding to the inelastic processes (oc F(eV±2J)). 
This is the consequence of the presence of the two-particle 
cigenstates 1 2 1 1 — 1) and |2110) which allow the quantum 
dot to make inelastic transitions using different paths. As 
usual it leads to the dependence of the probabilities of 
such transitions on phases of the tunneling amplitudes. 
This fact differs problem under consideration from the 
standard one [l7j in which only elastic terms depend on 
phases of the tunneling amplitudes. Low-temperature 
current-voltage characteristics are non-linear but they 
became linear if the gap |A — 2J\ tends to zero. One 
can utilize this fact to observe the quantum dot ground 
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state transition point experimentally. 

To realize effectively two-level quantum dot one can use 
any system with doubly degenerate levels and study it at 
voltages and temperatures much less than level spacing. 
For example, it can be a carbon nanotube with orbital 
degeneracy of levels [2£| or a two-dimensional electron 
gas in Si(001)-MOSFET [H and Si/SiGe heterostruc- 
tures [33l | . Also one will observe two-level quantum dot 
behaviour of inelastic co-tunneling current for arbitrary 
quantum dot in the regime |V|, T ~ | J — Jg\ « J, A. 

The authors thank A. Ioselevich for useful discus- 
sions. The research was funded in part by the Russian 
Ministry of Education and Science under Contract No. 
P926, RFBR Grants No. 09-02-92474-MHKC, the Coun- 
cil for grants of the Russian President Grant No. MK- 
296.2011.2, the Dynasty foundation and RAS Programs 
"Quantum Physics of Condensed Matter" and "Funda- 
mentals of nanotechnology and nanomaterials" . 



Hamiltonian (|A1[) commutes with the total number of 
electrons N, the total spin square S 2 , and S z . Therefore, 
it is convenient to work in the basis of two-particle states 
\NSnS z ) presented in Fig.Q] Then Hq D can be written 
as a 16 x 16 matrix. All the states except those with 
N = 2, S = 0, S z = and with N = 3, S = 1/2, 
S z = ±1/2 are eigenstates of Hamiltonian (|A1[) . Their 
energies are 

-E'oooo = 0, 

= #l|l-§ = £ 1> ^liQi = ^liO-i = £2, 
E^lll = E2H-I = -EimO = e l + e 2 + U\ 22 i — UxX22, 

-E4020 = 2ei + 2e2 + Uxxxx + U 2222 + 4J7i22i — ^U\x 22 - 

(A5) 

The states |2010), |2020), and |2000) are mixed and 
Hamiltonian (|Alj) projected onto these states is given 
as 



Appendix A: General form of Hamiltonian for a 
two-level quantum dot 

Although there is vast body of studies on spectra in 
a few electron quantum dots [34 |. in this appendix we 
present the general expression for Hamiltonian for a two- 
level quantum dot and discuss under which assumptions 
it can be simplified to Eq. ([3]) . We start from the follow- 
ing Hamiltonian 



HQD = 2_j eat ^a a dacr + Hint 



(Al) 



where a = 1,2 denotes the orbital single-particle levels 
and 



1 



a 2 a 3 Q 4 rf ai (jj rf 02 ff2 rfa 3 ff 2 ^Qj 1 (A2) 



is the interaction part of Hamiltonian. The matrix ele- 
ments of interaction arc defined as: 



U aia2a3ai = J drdr'(p* ai (r)ip* a2 (r')U(r - r') 

x (p a3 (r')ip a4 (r) (A3) 

where U(r) = e 2 /(er) is the Coulomb potential. Pro- 
vided the time-reversal invariance is preserved only six 
matrix elements U, 



aiQ 2 a3a4 



are independent: 



UXXXX: 








UXXX2 


= fll21 


— Ux2XX 


— l/2111 


f 1122 


= C^2211 


= ^1212 


= t/2121 


"l221 


= C^2112: 






(7i222 


= &2122 


= C^2212 


= (72221 


^2222- 









(A4) 

Here matrix elements f/im, U 2222 an d t/1221 describe 
direct Coulomb interaction whereas Uxxii corresponds to 
exchange energy. 



Hx = (ei+e 2 )l + Vi 
where 1 denotes the unit matrix and 



V x = 



U\ 22 x + Uxi22 
-V2UXXX2 
-v2Ux222 



-V2Uxxx2 -V2U! 222 

-A + dilii UxX22 
U\X22 A + C/2222/ 



(A6) 



(A7) 



The states |3i2|) and |3^1^) are also mixed and Hamil- 
tonian (| Al|) projected onto these states can be written 

as 



H 2 = {ex + e 2 + 2?7i22i - U 1122 )l + V 2 



where 
V 2 



ex + Uxxxx 
—Uixi 2 — U\ 222 



-UxXX2 — C^1222 \ 
e 2 + 1^2222 / 



(A8) 



(A9) 



Hamiltonian (|A8|) describes also mixing of the states 
|3±2-|) and |3§1 - §>. 

As an example, let us consider a quantum dot 
fabricated in two-dimensional electron gas in Si(001)- 
MOSFET structure [HJ. In such quantum dots electrons 
can occupy states in two valleys which remain from six- 
fold degeneracy of bulk Si. Assuming the level spacing 
due to spatial confinement to be large as compared with 
the valley splitting we have only two low-energy orbital 
states (symmetric and anti-symmetric): 

Qz 

<Px(r) = V2cos —(p (z)(pj_(p), 
Qz 

tp 2 (r) = V2sm—(p (z)ipj_(p) (A10) 

where z denotes the coordinate perpendicular to the 
plane in which two-dimensional electron gas is formed, p 
the in-plane coordinate vector, and r = p + ze z . The 
vector Q = (0,0, Q) corresponds to the shortest dis- 
tance between the valley minima in the reciprocal space: 
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Q oc 1/aiat, with ai a t being the lattice constant. Wc 
choose the envelope function 4>o(z) to be normalized and 
symmetric. We assume that the ground state eigcnfunc- 
tion tp±(p) of the confinement potential which creates a 
quantum dot is also normalized and symmetric. Then, 
as one can check the matrix elements Um 2 = ^1222 = 0. 
Also wc assume that the following conditions 

Q- 1 <£l z <s:l±_ (All) 

are satisfied. Here l z and l± are typical lengths associated 
with the functions ifio(z) and ip±(p): I" 1 ~ J dz (p%{z), 
I J 2 ~ J dpip^(p). Then neglecting exponentially small 
in Ql z contributions of a type J dz (pf)(z) cos(2Q z) , we 
obtain 

Chilli = £^2222 = U + J, C/l221 = U — J, U\\22 = J 

(A12) 

where 

U= l -J drdr'U{r-r') V l{z)vl(pWl{z')vl(p') 
f/ii22 =\U{Q) J dziptiz) J dpcpi(p). (A13) 

Here U(Q) — 4ire 2 /(eQ 2 ) stands for the Fourier trans- 
form of the interaction potential. The direct Coulomb 
energy can be estimated as U ~ e 2 / (el±) and is just the 
charging energy E c . The exchange energy Un 22 can be 
estimated as Un 22 ~ e 2 /{Q 2 l z l\) ~ U/(Q 2 IJ±) < U. 

Therefore, the states with two electrons on the quan- 
tum dot can be described by the following Hamiltonian: 

H 3 = ]T e a d) a<J d aa - JS 2 + A c TtT (A14) 

oca 

where T = ^ d a ^d a ± and A c = J. The last term in 
Eq. (|A14j) describes superconducting correlations. How- 
ever, since the interaction in the Cooper channel is repul- 
sive (A c > 0) one can expect that it will rcnormalize to 
zero due to virtual transitions to high energy levels Q. 
Then Hamiltonian (|A7j) coincides with Hamiltonian ((3]) 
projected to the states with two-electrons on the quan- 
tum dot. 



Appendix B: Explicit expressions for x' s 

We present explicit expressions for the x' s which one 
needs to know in order to calculate the co-tunneling cur- 
rent (fT5j) . Also we present expressions for the averages in 
Eq. ([SJ in terms of the exact two-particle correlators for 
the isolated dot (Hqjj) and the Green functions of elec- 
trons in the leads which are used to calculate x' s - Wc 



obtain 

(X*(t)H T (ti)H T (h)H T (t 3 )) = 

(da lt da 2 tidl :st2 d a4t ^) T 1 4s 2 G l ^ i ^ 4 (t,t 3 )G r ^ 2 (t2,ti) 

- (da lt d a2 t 1 dl t3t2 d a4 t 3 ) Tl23iG l ^p 2 (t, h)Gp>p 4 (t 2 ,t 3 ) 

+ (d^tdi^da^da^ts) n 32 4G l ^ 1 p s (t,t2)Gp 2 > p A (t 1 ,t 3 ) 

- (d^tdi^da^daita) T 1 4 23 G l ^ i p 4 (t,t 3 )G r p 2 > p 3 (t 1 ,t 2 ) 

+ (di 1 tda2t 1 d a3 t 2 dl ilt3 ) T 12i3 G l ^ il32 (t, h)Gp2 h (t 3 , t 2 ) 

- (d' ! ai tda2t 1 d a3 t 2 dl lit3 ) n 3 42G l ! j > i p 3 (t,t2)G r p<p 2 (t 3 ,ti), 

(Bl) 

where (...) = Tr . . . e"^ H « D / Tr e -^ D , a k = {a, a}, 
njkl =tp iai tl f3 ja ?hcc k tr f) im and d a t = d a (t). Next 

(da 1 t 2 dl 2t d a3tl d'l it3 )T 23 4 : iG l ^' 2l33 {t,ti)G r p<p 1 (t3,t2) 

- (da 1 t 2 da 2t d a3 t 1 d'l it3 )T 21 4 3 G l l ^ j3i (t, t2)G r p<p 3 {t 3 , ti) 

+ (da^da^dast^a^T^lsG^^it^^G^p^^x) 

- (da^da^dast^a^T^liG^^it^^G^^itx,^) 

- (da 1 t 2 dl l2t dl l3ti d a4 t 3 )T 2 4 31 Gp > 2 p 4 (t, t 3 )Gp 3/3i (tl,t 2 ) 

+ (^Qit2^ 2 t d a3tl^ Q 4t3) T 2134G^ (3i {t,t2)G r ^p i (t 1 ,t 3 ), 

(B2) 

{HT{t3)H T {t2)X\t)H T (ti)) = 

(d ai t 3 d a2 t 2 dl l3t d' ! aitl )T 31 4 2 Gp < 3j3i (t, t 3 )Gp 4/32 (ti,t 2 ) 

- (d ait3 d a2t2 d ast d aitl ) T 32 4i G 1 ^^ (t, t 2 )G r p2 /3i (h,t 3 ) 

~ (^ 1 t3 rf a2t2^i 3 t^a4tl) T 3412G^ /34 (t, tljG^^ (t 3 , t 2 ) 
+ (da 1 t 3 da 2 t 2 dl 3t d a4 t 1 )T 23 14 : Gp > 2 p 3 (t, tl)Gfap A (t3, t 2 ) 

- (d ait3 d a2t2 d a3t d aitl )T 3 i 24 Gp < 3/3i (*,*3)G^ 4 (<2,<i) 

+ (d ai t 3 da 2 t 2 dl 3t d a4tl )t 34 21 Gfa0 4 (*j 1 1 )G f fap 1 (h ,t 3 ), 

(B3) 

(fr T (ti)A-t(t)flT(t2)flr(t3)> = 

{d ai t! da 2t d a3 t 2 d ait3 )r 2 34i G%p 3 (t, h )G , ^ (3l (t 3 ,h) 

- (d ai t! da 2t d a3t2 d ait3 )t 2 i4 3 G^ /3i {t,t\ )Gjg^g 3 (*3 ) h ) 
+ (d^^da^da^da^T^lsG^pJt^^G^pJtx,^) 
~ (da^rda^das^da^T^uG^pJt.t^G^pJti, t 3 ) 

- {d ai t 1 d\ 2t d\ 3t2 d ai t 3 )T 2 4 3 iG l p > 2 p i (t, t 3 )G r p 3/3i {t 2 ,ti) 

+ (d ai t 1 d'l 2t d'l 3t2 d ai t 3 )T 213 4G l l ^ l3i (t,ti)Gp^p t (t2, t 3 ), 

(B4) 
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(HT^HTit^X^^HTih))) = 

(d ai t 3 d a2tl d a3t d ait2 )T3U2G l ^p i (t, h)G~p2p 2 (t 2 , ti) 

da2t 1 d a3t d a4 t 2 ) T 3412G fcpjt, ^2)G^ 2 (t3, t\) 
da2t 1 d a3t d ai t 2 ) T 3214G fop 2 (t, tl)Gp^p i (t3,t2) 

- (d ai t 3 da 2tl dl 3t d ai t 2 )T3i24G l ^ i p i (t, ts)G^p 4 
+ (da 1 t 1 da 2t d a3 t 2 dl it3 )T234iG l l j > 2 p ;i (t, t 2 )G r ^p i (t 3 , t\) 
+ (d ai t 3 da 2tl dl 3t d ai t 2 )T 3 42iG l ^ 3l34 (t, t 2 )G r l ^ l3i {ti,t 3 ), 

(B5) 



(fT r (ta)i^ti)*t(t).ffr(t3)) = 

(d ai t 2 d a2tl d ol3t d ait3 )T3U2G l ^p i (i, < 2 )G^ / 3 2 (*3, *i) 

- (d ai t 2 d a2tl d a3t d ait3 )T324iG l £ s p 2 (t, t 1 )G r p2 l3i (t 3 , t 2 ) 

da 2 t l d^ a3t d ai t 3 ) T 34l2G l p 3l}i {t 1 t3)G r l ^ f}2 {t2 1 ti) 

+ «* 2 

da 2 t 1 d a3t d ai t 3 )T3214G ^p 2 (t, 1 1 )G^ 4 (t 2 , ^3) 

- (d ait2 d a2tl d a3t d ait3 )T3i2AG l p 3 p 1 (t, t 2 )Gp 2 > p 4 , t 3 ) 
+ (da 1 t 2 da 2tl dl 3t d ai t 3 )T342iG l p 3 p 4 (t, t 3 )G r p 2l3i {h, t 2 ), 

(B6) 

(H T (t3)H T {t2)H T (h)X\t)) = 

(d ai t 3 d a2 t 2 d u 3 t 1 d u i t) T 423iG p ij}2 (t, t 2 )G r l ^p i (ii, i 3 ) 

- (d ai t 3 rfa 2 t 2 d a3tl d Q , 44 )r4i3 2 G^ /9l (t, t 3 )Gp 3/32 < 2 ) 

- (d Ql t 3 di 2t2 da3t 1 di 4t )T43 2 iG'^ (33 (t,ti)G'J 32 ; / g i (< 2 ,i3) 
+ (rfQit3 rf i 2 t 2 rf a 3 ti^ 4t )r4i 23 G'^ (3i (t,i3)G'J 32 > (33 (t 2 ,<i) 

+ (da 1 t 3 da 2 t 2 d a3 t 1 dl l4t )T43 1 2G l ^p 3 (t,t 1 )G r p> l32 (t3,t2) 

- ( d a!t 3 da 2 t 2 d a3tl <4 4 1) 7"4213 Gjg^ (t, t 2 (t 3 , 1 1 ) . 

(B7) 



(da 1 t3d a2 tda 3 t 1 d ait2 )T234lGfa j33 (t,t 2 )Gp<p 1 (t3,ti) 

- (rfait3 rf l 2 t rf a3tl^i 4 t2) T214 3G^ / 3 1 (i,i3)G^ /33 (t3,ti) 

d a2t da 3 t 1 d ai t 2 )T2413G ^ 2j34 (t, £ 2 )G^ 3 (£3 , ti ) 

- ( rf l 1 t 3 ^i 2 t^a3tl^a4t2) T 2314G^ /33 (t,ti)G^ /34 (t 3 ,t 2 ) 

- (d ai t 3 d'l 2t d'l :jti d ai t 2 )T243iG l ^ 2l34 {t, < 2 )G^ 3 c /3i (ii, t 3 ) 
+ (da 1 t 3 dl l2t d'l :jti d ai t 2 )T2i34G l l ^ l3i (t, i 3 )G^ 4 * 2 ), 

(B8) 



The non-interference contributions to %'s are given as 



X<'> = A 2 Re f dh f dt 2 f dt 3 ( 
J —00 J —00 J —00 

+ 9a 1 9a 3 (da 1 t 3 dl ll tda 3 t 1 d a3 t 2 ) 

+ 9a 1 9a 3 (da 1 t 2 dl ll tdl l3 t 1 da 3 t 3 ) 

~ 9a 1 9a 2 {da 1 t 3 d\ X2t2 d\ xlt d a2 t 1 ) 

+ 9a 1 9a 3 {da 1 t 1 d\ xit d\ 13t2 d a3 t 3 ) 

+ 9a 2 9a 1 (da 1 t 3 da 2 t 1 dl l2t d ai t 2 ) 

+ 9a 3 9a 1 {da 1 t 3 d\ xltl d\ X3t d a3 t 2 ) 

+ 5i 2 ffa 1 (^ 1 t 2 ^a2tl^l 2 t^ait 3 ) 

— 9a 1 9a 2 {da 1 t 2 d\ X2tl d\ xlt d a2 t 3 ) 

+ 9a 1 9a 2 {da 1 t 3 d\ X2t2 d a2 t 1 d\ xlt ) 

+ 9a 3 9a 1 (da 1 t 3 da 1 t 2 d a3 t 1 dl l3 t) 

~ gi^aMaitJaztzda^dirf)) , (B9) 



X > ' < =A 2 Re / rftj / di 2 / dt 3 ( 
J —00 J —00 J —00 

~^9 ai 9a 2 (da 1 td a2 t 1 dl l2t2 d ait3 ) 

~ >r 9a 1 9a 3 (da 1 tda 1 t 1 d a3 t 2 d a3 ta ) 
_ 3i 1 5Q 2 ( d l 1 t rf a 2 tirfait 2 rfi lt3 ) 

+3a 2 5ai (da 1 t 3 dl l2t d a2 t 1 dl ll t 2 ) 

~9a 2 9 ai (da 1 t 3 dl l2 tda 1 t 1 da 2 t 2 ) 

+9 l a2 9u 1 {du 1 t 2 d { a2t d a2 t 1 d { ait3 ) 

~9a 2 9 ai (da 1 t 2 dl l2 tda 1 t 1 da 2 t 3 ) 
~ >r 9a 3 9a 1 (da 1 t 3 dl 11 t 2 dl 13 tda 3 t 1 ) 

+9 l a2 9a 1 {d ai t 1 d { a2t d a2 t 2 d { ait3 ) 
~9a 2 9ai (dai_tid a2t d ait2 dot 2 t 3 ) 
+3a 3 3ai (da 1 t 3 dl lltl dl >l3t d a3 t 2 ) 

+9 l a3 9 r ai ( d 

ait 2 da ltl C?jj 3t C?Q, 3 t 3 ) , (BIO) 
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A Re 



dti 



to 



(It 



{ 



+9 l ai 9 r a2 

I r 
9ol\ 9a.2 

+g l a2 g r ai 

I r 
9a.2 9ai 

+9 l a2 9 r ai 

I r 
9a-2 9ol\ 

I r 
9 at j, 9a 1 

+9 l a2 9ai 

+9 l a2 9 r ai 
-9 l a2 9 r ai 

l T 
9a 3 9oi\ 



daitda^! dL 3 t 2 da 3 t 3 
^ait^a 2 ti < ^a 1 t 2 d a 
Uj a 1 t u a 2 t 1 Ll a 2 t 2 «cti 



2t 3 / 
t,) 



d~,„±(lrvi i-i dr\ 



u a 1 t 3 a a 2 t u ait 1 U a2 t 2 

jt jt J J 
a a 1 t 3 u a 2 t u a2*l u ai*2 

jt jt J J 

U Qlt 2 U, a 2 t U '"l t l U ' tt 2t3 



jt jt J J 
u ait 2 u a2t a 2*l u aii 3 



'iit3 a ai*2 a a 3 t a a3ti 



da 1 t 3 d a2 t 2 dl X2t d ai t 1 ) 

da 1 t 1 da 2 tda 1 t 2 d a2 t 3 ) 
da 1 t 1 da 2 tda 2 t 2 d ai t 3 ) 

da 1 t 3 d ai t 1 dl 13t da 3 t 2 ) 



-9a 3 9a 1 (da 1 t 2 da 1 t 1 d' ! a3t d 



(Bll) 



r<'< ee A 2 Rc 



(its 



'2 



9 cui 9 as 

+9 l ai 9 r a2 

I r 

9a.-2 9ol\ 

I r 

9 011 9a>3 

+9 l ai 9 r a2 

I r 

9012 9 Oil 

+9 l ai 9 r a2 

l r 
~9a 2 9 ai 

+9 l a2 9 r ai 

l r 
9a±9a 2 



da 1 t 3 d] }llt d a3 t 1 d a3t2 ) 
da l t 2 d\ xit d a3 t 1 d ast3 ) 

da 1 t 3 d a2 t 2 dl ll tdl l2 t 1 ) 
da 1 t 3 d a2 t 2 dl l2t dl lltl ) 

daitid ait da 3 t 2 d a3 t 3 ) 

d ai t 3 d a2 t-L d ai t d a2t2 ) 
d ai t 3 d a2 tid a2t d ait2 ) 
da 1 t 2 d a2 t 1 d^ ollt d rX2t3 ) 
d ai t 2 d a2 t 1 d a2t d ait3 ) 



d(x^t 3 d. 



a 2 t 2 dl ll t 1 d a2 t) 
ta 1 t 3 d a2 t 2 dl l2 t 1 d jf ai t) 



-gl 3 9 r ai (d ai t 3 dt lt2 d a3tl di 3t }) . (B12) 
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